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ABSTRACT 

New molecular (^^CO J=3— 2) and dust continuum (450 /um and 850 /um) maps of the 
NGC7129 star forming region are presented. The maps include the Herbig Ae/Be star LkHa 
234, the far-infrared source NGC 7129 FIRS2 and several other pre-stellar sources embedded 
within the molecular ridge. The data are complemented with C^^O J=3— 2 spectra at several 
positions within the mapped region. Both the submillimeter and the ^"'CO emission show a 
similar morphology, displaying a sharp boundary towards the cavity. The submillimitcr maps 
also reveal a second source, SMM2, which is not clearly seen in any earlier data set. This is 
either a pre-stellar core or possibly a protostar. Also, the highest continuum peak emission 
is identified with the deeply embedded source IRS 6 a few arcseconds away from LkHa 234. 
These new 450 and 850 /xm observations are combined with previous continuum observations of 
the three compact far-infrared sources in the field, in order to make fits to the spectral energy 
distributions and to obtain the source sizes, dust temperatures, luminosities, and masses. For 
nine positions where both ^^CO and C^^O spectra are available, gas masses have been obtained 
and compared with masses derived from the continuum fluxes. The masses are found to be 
consistent, implying little or no CO depletion onto grains. The dust emissivity index is found to 
be low towards the dense compact sources, /3 ~1 — 1.6, and high, (3 ~ 2.0, in the surrounding 
cloud. 

Subject headings: ISM: dust, molecules - submillimeter - stars: formation, pre-main sequence - 
stars: individual: LkHa 234 
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Introduction 



NGC 7129 is a reflection nebula seen against a molecular cloud and estimated by Shevchenko et 



al. (1989) to be at a distance of 1.25 kpc. There are many signs in this region which support the idea 
of triggered and possibly ongoing star formation: the association of a very young cluster formed by the 
B3-type stars BD +65 1637 and BD +65 1638, and the B5/7-type star LkHa 234, several embedded infrared 



sources (Harvey et al. 1984; Weintraub et al. 1994; Cabrit et al. 1997) as well as reflection nebulae and 



Herbig-Haro (HH) objects (Hartigan & Lada 1985; Miranda et al. 1994). There are also several molecular 



outflows in the region: one associated with an optical jet near LkHa 234 (Ray et al, 1990; Edwards & Snell 



1983| ), another bipolar outflow near the far-infrared source NGC7129 FIRS 2 and one may be driven by the 



T Tauri star V350 Cep (Hartigan & Lada 1985; Goodrich 1986; Miranda et al. 1994) 



A molecular cavity is revealed in both molecular and infrared maps (Bechis et al. 1978; Bertout 1987), 
but its origin is still unknown. The location of the two bright stars, BD +65 1637 and BD +65 1638, inside 
the cavity, suggests three mechanisms by which the molecular cavity could have been created: (i) through 
the advancement of a photon-dominated region (PDR) around those stars, in which the molecular gas was 
dissociated by their intense UV radiation, (ii) by an expanding shell associated with stellar mass loss, or 



(in) by radiation pressure acting on grains, together with gas-grain drag (Bechis et al. 1978). The two 
mid-B type stars in the cavity, BD +65 1637 and BD +65 1638, are also the oldest in the cluster. Therefore 
stellar mass loss may have been very effective at earlier times and the winds from those two stars may have 
triggered the formation of other bright stars in NGG 7129, including LkHa 234 and SVS 13. As for what 
mechanism is dominating at present time, the question remains open. PDR models can marginally explain 
the line fluxes of [01] and [CII] emission in the region, although the C-type shock models cannot be ruled 
out (Lorenzetti et al. 1999| ; Giannini et al. 1999). 



This paper employs new, high resolution submillimeter observations of the molecular ridge in dust 
continuum emission, ^^CO line emission, and C^^O line emission. One goal of this study is to improve 
the present picture of the structure of the ridge and the ridge-cavity boundary. Another goal is to obtain 
the physical properties of the compact far-infrared sources in the ridge. A third aim is to obtain the 
gas/dust ratio, thus determining whether CO isotopomers are a good surrogate for molecular hydrogen, 
or whether CO has substantially frozen out onto dust grains. The final goal is to search for variations in 
the dust opacity with position through the field. The mapped region includes the dense molecular cloud, 
a photodissociation region, and a powerful outflow. We might expect processes such as grain coagulation 
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and mantle evaporation to have a strong spatial dependency here and for this to be reflected in spatial 
variations of the dust emissivity index, (3. 



2. Observations and data reduction 

2.1. ^^CO and C^^O observations 

Observations of ^'^CO and C^^O J=3— 2 were made with the James Clerk Maxwell Telescope, on 
Manna Kea, Hawaii, in August 1994 and November 1997 respectively, with the receiver B3i and the Digital 
Autocorrelation Spectrometer (DAS). The observations were made by position switching, with the reference 
position 10' East of LkHa 234, which was determined to be free of significant CO emission. The main 
beam efficiency (?7{,m) was 0.58 and the half-power beam (HPBW) was ~ 14" for both ^^CO and C^^O 
observations. 

A grid map with the spacing of 7'.'5 x 7'.'5 was taken in the ^^CO J=3— 2 (330.58 GHz) rotational 
transition. The map was centered near the star LkHa 234|, at coordinates a(2000) = 21'' 43™ 06". 08 and 
<5(2000) = 66° 06' 56'.'29, and extending from -100" to +100" in a and from -150" to 75" in 5 offsets. 
The instrumental velocity resolution was 0.28 km s^^. C^^O J=3— 2 (329.33 GHz) spectra were taken at 
nine positions around the star LkHa 234 and near the molecular ridge. Figure 1 shows the ^^CO integrated 



intensity (/ TbcIv) contour map superimposed on a K' image of the NGC 7129 region from Hodapp (1994) 
The nine C^^O offsets with respect to the star LkHa 234 are shown with cross symbols. 



2.2. Submillimeter continuum observations 

Two deep 850 /im and 450 /xm-maps of LkHa 234 and the bright-rimmed cloud ridge southwest of 
LkHa 234 have been obtained using the Submillimeter Common User Bolometer Array, SCUBA (Holland 

et 



1999), on the 15m James Clerk Maxwell Telescope, on Mauna Kea, Hawaii, on October 15, 1997. 
The weather conditions were dry and stable with an 850 /im zenith optical depth of ^ 0.19. The two 
overlapping fields were obtained in jiggle-map mode with a chop-throw of 150" in Right Ascension with a 
total integration time of 85 minutes/field. For a description of SCUBA and its observing modes, see Holland 



^The astrometric position for LkHa 234 given by Clements & Argyle (|l984|) is a(2000) = 21'' 43" 06^802 
and 5(2000) = 66° 06' 54"5. 
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et al. (199£). Unfortunately the maps were taken without usable calibration and pointing observations 
and we have therefore reduced the data by determining the HPBW of the telescope from observations of 
Uranus from other nights in September and October with similar sky conditions. The same nights were 
also used to determine gain conversion factors for the 850 yum and 450 /im filters. The HPBW was found 
to be ~ 15'.'6 X 13'.'5 at 850 fim and 9'.'1 x 7'.'8 at 450 fim with the beam broadened in the chop direction. 
The submillimeter position of LkHa234 (Table 3), which has an accuracy of 1", was taken from three 
short-integration maps interlaced with pointing observations in SCUBA commissioning time during Spring 
1997. These maps have not been added in the final data set, because they were taken with a 120" Azimuth 
chop and there is clear evidence that the array was chopping onto emission. However, if we compare the 
total flux density of LkHa 234 (i.e. background subtracted flux density) it agrees very well with data taken 
in October. 

We also include a map of the H2O maser and bipolar outflow source NGC 7129 FIRS2, which was used 
as a test source during SCUBA commissioning time (SandeU |T997| ). A composite map at 850 /xm, including 
all the three overlapping fields is shown in Figure 2. 



The basic data reduction was done with SURF (Jenness & Lightfoot 199E). Even with a 150"-chop it 
is more than likely that we may have chopped onto some emission in both fields, but it does not appear to 
affect the morphology of the strong dust ridge in which LkHa 234 is embedded. Some negative emission is 
detected west of LkHa 234, which confirms that there is some emission in the off-position. It is therefore 
likely that emission east of the ridge is also affected, which is why the 850 and 450 /xm maps appear 
somewhat different. To map a cloud ridge in beam switch mode is not very advisable, but at the time 
SCUBA was only available in jiggle map mode. Although chopping into the cloud cannot be avoided, as 
long as the off-source emission is faint and relatively uniform, the morphology of the cloud ridge should not 
be affected. However, since the beam size and error lobe contribution differ between 850 /im and 450 /xm, 
this leads to systematic errors in the integrated intensities (Section 3.2). 

The peak fluxes are uncertain to about 10% for 850 /im and 20% for 450 /xm. The final pointing 
corrected and calibrated maps have an rms noise of ^ 15 mJy/beam and ^ 120 mJy/bcam for 850 and 
450 /im, respectively. These maps were written out as FITS files and read into the MIRIAD package (Sault, 



Teuben, & Wright 1995) for further analysis 



3. Analysis 
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3.1. i^co and C^^O Emission 

A sharply bounded molecular cavity is seen in ^'^CO 3—3—2 emission westward of LkHa 234, most 
probably excavated by the B3-type star BD +65 1637. The lack of molecular emission and the ridge 



morphology were also noted in the CO J=l— and NH3 observations (Bechis et al. 197S; Fuente et 



1998; Glisten & Marcaide 1986). However, our high resolution (half-beam sampling with a 14" beamsize) 



observations resolve the ridge structure in much more detail. 

There are four ^^CO peaks along the ridge, visible in Figure 1 as contour peaks. These ^'^CO structures 
are assumed to be physical clumps and named NGC7129 COl, NGC 7129 C02, NGC 7129 C03, and NGC 
7129 C04. The positions of these four ^■^CO peaks, together with their peak antenna temperatures, line 
widths, and line center velocities for the ^■^CO emission are given in Table 1. NGC7129 COl coincides 
within errors with LkHa 234 (position 5 in Table 2), NGC 7129 C04 with position 7 and the other two 
clumps, NGC7129 C02 and NGC7129 C03 are in located close to the positions 6 and 3, respectively. 

Perhaps the most striking feature of the ^'^CO map is the sharpness of the transition between the 
molecular cavity to the west and the molecular cloud to the east. ^^CO line strengths drop by factors of 20 
or more across the boundary. The narrowness of the transition region between atomic and molecular gas has 
implications for the three dimensional topology of the region. There seem to be three possibilities: Firstly, 
the boundary region may be a sheet which happens to be perpendicular to the line-of-sight. Secondly, the 
molecular cavity may be a bowl, which is surrounded by molecular gas. Thirdly, the molecular cloud may 
take the form of a rather narrow ridge, which is partially surrounded by atomic gas. The first possibility 
is improbable, while the second should lead to a broader transition region (molecular to atomic). We favor 
the third picture, in which the dense portion of the molecular cloud consists of a clumpy ridge with a line-of 
sight extent comparable to its extent in the plane of the sky (roughly 0.5 pc). 

Another remarkable feature which can be seen in Figure 1 is the narrow filament of emission seen 
in the K' image (Hodapp 1994). It extends roughly south from LkHa 234 and then curves to the west, 
keeping a nearly constant surface brightness. In an unpublished observations using an IR camera with a 
circularly variable filter on the CFHT, Mitchell and Nadeau found this filament to be also a source of strong 
emission in the S(l) v= 1—0 line of II2. The optical filament follows closely the outermost ^'^CO contours, 
which separate the cavity from the molecular cloud. Can this close vicinity be explained by standard PDR 
models? The presence of the photodissociation region is a natural explanation for both the sharp molecular 
boundary and the line-emitting filament. The H2 S(l) v== 1—0 line emission is then due to fluorescence. 
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The PDR is very likely produced by UV radiation from the early B star BD+65 1638. However, the 
detection of all three types of emission: optical, CO and H2 v= 1—0, at the same location, is in contrast 



with predictions of steady state stationary PDR models (HoUenbach & Tielens 1997 ). These models predict 
a more gradual succession of emission layers as one moves away from the UV ionization front, and where 
the H2 v~ 1—0 transition layer would be separated by about 10" from the molecular CO emission. On the 
other hand, the observed morphology may be explained in the framework of nonstationary PDR models 



(Bcrtoldi & Draine 1996), which show that the inhomogeneities in the molecular cloud surface can lead to 



the merging of the ionization and the dissociation front. 

Because ^^CO is likely to have an appreciable optical depth, C^^O 3—2 spectra were obtained at several 
positions. The positions for which we have C^^O spectra are marked by the cross symbols in Figure 1. A 
comparison between the ^^CO and C^^O spectra at the nine offset positions is shown in Figure 3. Positions 
are given as offsets in arcseconds from LkHa 234. The line intensities are given in antenna temperatures, 
Ta*, after corrections for losses due to the Earth atmosphere and the telescope. For all spectra only the 



linear baselines were removed, using the SPECX software package (Padman 1992) 



Near LkHa 234, the lines are broadened due to contributions from the outflow. Position 1 (Table 2) is 
near the peak of CO emission in the redshifted outflow lobe. Both the ^"^CO and C^^O spectra are double 
peaked at this position. The value of the line ratio shows that C^^O is optically thin. The two peaks 
therefore represent two real kinematic components rather than a single line suffering self absorption. The 
peak at —10 km s''^ is the ambient gas at this position, while the peak at —7 km s^^ is the postshock gas. 

Several line parameters have been calculated at each offset, assuming local thermodynaniical 
equilibrium (LTE) conditions. A value of 10 was used for the abundance ratio Nisqq/ Nqi8q, estimated 
from the partial ratios Nci8o/Nh2 2 x 10^^ (Frerking, Langer & Wilson 1982| ) and the standard value 



Ni-iQol Ni-ico ^ 56 (Wilson and Rood, |1994D . The intensity ratio, ^^CO/C^®0, was used to find the optical 
depth, assuming a common for the two transitions. The value of Tg^, was computed from the LTE 
expression relating line strength (given by the radiation temperature), line optical depth r, the excitation 
and radiation temperatures: 



T =^ I (1) 

k \n[l + {l-e-^)hv/kTR]' ^ ' 

In the above relation, the radiation temperature, Tr, is related to the antenna temperature, T* a, by: 
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Tr = ^, (2) 

VbmJ 

where the beam efhciency 7yf„„ was 0.58 and the fiUing factor f was assumed to be 1. We will further assume 
that the gas kinetic temperature is equal to the excitation temperature, but caution that a density greater 
than 10"* cm""^ is required for this to be true. The results of the LTE analysis are shown in Table 2. For 
LkHa 234, the excitation temperature is found to be about 30 K, in agreement with previous values of 20 



— 30 K obtained from lower resolution (70") CO observations (Bechis et al. 1978) and with the 20 — 25 K 
obtained from NH3 (40") observations (Giisten & Marcaide [l98^ ). 



At position 6 (Table 2), the gas has an excitation temperature of ~ 34 K, which is higher than the 
typical value for a dark core, but lower than the previous value of 63 K obtained by Mitchell and Matthews 



( 1994 ) from ^^CO and ^^CO data. These two temperatures can be reconciled, taking into account that the 
isotopomcrs from which they have been derived (^^CO for 63K and ^■^CO for 34 K) have different optical 
depths and hence probe different regions of the gas. With its lower optical depth, ^■^CO will probe deeper 
into the PDR, where the temperature is lower. 

The three last entries in Table 2 represent positions increasing in distance from the cloud boundary. 
The excitation temperature is found to decrease with distance from the boundary, from 38 K through 31 
K to 21 K. This decrease is expected as the stellar radiation is increasingly attenuated with distance into 
the molecular cloud. At position 1 (Table 2), which is farthest from the cloud boundary, the excitation 
temperature of the ambient gas is lowest, ~ 13 K. 

Since the clumps have approximately the same size as our beam and since the beam does not always 
include a clump of gas, it is reasonable to calculate the masses of the gas within the beam (14"). The 
results, including the fractional He abundance and assuming a distance of 1.25 kpc to the source, are shown 
in Table 2. The NGC 7129 COl clump which coincides with LkHa 234, contains 10.7 Mq. The masses 
within the beam, at the two offsets adjacent to the clump matching the pre-stellar source LkHa 234 SMM 2 
(positions 3 and 6 in Table 2) are 12.8 and 7.2 Mq, respectively. However, those values should be regarded 
as upper limits since the temperatures are probably lower, as suggested by the more optically thin NH3 
data. 

These results also depend on the accuracy in the distance determination. The distance estimate of Ikpc 
to NGC7129, reported by Racine ( 1968| ), has been most commonly used. In our analysis we have adopted 



the value of 1.25 kpc, determined by Shevchenko et al. ( |1989| ) on an improved photometry and spectral 
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classification of 35 stars. However, if the PDR is associated witli to the Cephcus Bubble, as suggested by 
Abraham , Balazs & Kun ( 1999| ), then the distance to could be as low as 500—700 pc, and of course, the 
derived masses and luminosities would be significantly reduced. 

The masses obtained in this analysis will be discussed and compared to masses obtained from the dust 
emission, in Section 4 below. 



3.2. Dust Emission 

In order to compare the 850 and 450 /im dust emission with the molecular line observations, the 450 fj,m 
images need to be convolved to the same HPBW as the rest of the observations. However, since the JCMT 
telescope is far from perfect at 450 /im, one first needs to remove the error beam. A model beam was 
therefore constructed, by fitting three Gaussians to representative beam maps of Uranus. The main beam 
and the near error lobe were fitted with elliptical Gaussians, while the far error lobe was fitted with a 
circularly symmetric beam. The relative amplitude of the main beam and the near and far error lobes were 
found to be 0.925:0.07:0.005. The HPBW of the near error lobe was 37" x 26", with the major axis aligned 
in the chop direction similarly to the main beam (c.f. Section 2.2), and the far error lobe gave a HPBW ^ 
120". The deconvolution of images was tested with both MAXEN (MIRIAD's maximum entropy task) and 
CLEAN. Both worked reasonably well, but CLEAN is better in preserving flux, and was therefore adopted. 
We also constructed a model beam for 850 and and ran CLEAN on the 850 /im-images as well. The 
deconvolved images were restored with a symmetric Gaussian to a 14" beam. For 450 /im we also made 
maps restored to a HPBW of 8". Figure 4 shows the high resolution maps at 450 /im and 850 /im. 

Even CLEAN, however, is not ideal for recovering smooth extended emission, which could cause to 
underestimate the 450 /im emission more than at 850 /im, since the 450 /im-map is restored from an 8'.' 5 
resolution to 14", while the resolution in the 850/im-maps remains roughly the same. Emission in the off 
position will have an even more severe systematic effect. Since the error beam has a much higher amplitude 
at 450 /im than at 850 /im, any emission in the off position will affect the baselevel of the 450 /im-map 
much more than that at 850 /im. This leads to a systematic underestimate of the flux densities at 450 /im. 
This was tested by adding a constant flux level of 5% of the peak value at 450 /im to the "raw" 450 /im-map 
of LkHa234. The baselevel corrected map was run through CLEAN and restored to 14". When the flux 
densities in the 14" beam are compared for positions more than 40" away from LkHa 234, it is found that 
adding a zero-level of 5% of the peak flux density increases the flux level by more than a factor of two 
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compared to the uncorrected map. Since the base level can easily have an error of 5%, this means that 
the flux densities for faint extended emission are uncertain by more than a factor of two. In this case the 
450 fim flux densities are more than likely to be systematically underestimated. 

Two-dimensional Gaussian fits were then used to determine the size and total flux of the two "compact" 
sources in the field, i.e. LkHa234 SMM 1 and the protosteUar source, LkHa234 SMM2, - 28" North- West 
of LkHQ!234 SMMl. Here the fiux estimates are much more accurate, since the sources are compact and 
the emission from the surrounding cloud ridge can be substracted. LkHa 234 SMM 1 emission is found to 
be extended, with a size of ^ 5" (Table 3). This agrees well with earlier 800 /im and 450 /^m maps (Sandell 
& Weintraub 1994 ) as well as with 1.3 mm maps by Fuente et al. ( 199^ ) and Henning et al. (1998). The 



source LkHQ!234 SMM 2 is not obvious in the 1.3 mm data of Fuente et al. (1998) nor in the map published 



by Henning et al. (1998), but both 1.3mm maps show an extension towards SMM2 in good agreement 
with our SCUBA maps. By using the size derived from our SCUBA maps, we derive an integrated 1.3 mm 
flux of 0.2 ± 0.04 Jy for the map published by Fuente et al. ( |1998D . We are interested to derive the /3-index 
(where /3 reflects the change in dust emissivity or, equivalently, the change in optical depth, with frequency, 
K (X r oc v^^) and the total (gas -f dust) mass. For this, we have fltted the 1.3 mm data point and the two 
SCUBA data points with a simple isothermal model (see e.g. Sandell 2000| ). Since both IRAS and Kuiper 



Airborne Observatory (KAO) data (Bechis et al. 1978; Harvey, Wilking & Joy 1984; Di Francesco et al. 



199q) employ large beam sizes that are likely to include hot dust from the surrounding reflection nebulosity 



as well, we have not included FIR data in our least squares flt, although we constrained our LkHQ:234 



SMMl model to the peak dust temperature derived by Harvey, Wilking & Joy (1984). The results of the 
fit are given in Table 4 and also shown in Figure 5. LkHQ;234 SMM 2 is even less constrained. However, 
by looking at gas temperatures, either from our own CO observations or at the temperature estimates from 



the NH3 map by Giisten & Marcaide ( 1986 ), we find the gas temperatures to be < 30 K. Therefore the flt 
presented in Table 4 was constrained to 30 K. The mass we derive for the dust disk or envelope surrounding 



LkHa234 SMMl (Dent et al, |1989| ), is about half of that derived by Henning et al. ( |1998D after it is 
adjusted to the same distance. The difference is largely due to the difference in the adopted mass opacity 
for the dust. 

In order to compare the gas mass estimates derived from our ^^CO and C^^O data with dust mass 
estimates, submillimeter flux densities have been derived from our submillimeter maps convolved with a 
14" beam. The results are given in Table 5. If one uses the baseline un-adjusted 450 /im data and ignores 
the differences in dust properties and dust temperatures between the two embedded stars by assuming that 
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the dust is at the same temperature as the CO, the /3- values appear unplausibly low. If instead the baseline 
adjusted 450 /im data were to be used one would obtain a /3 of ^ 1.7. However, since we cannot accurately 
correct for the emission in the off position, we prefer to assume a /3-index and only use the flux density 
measurements at 850 /im, which are much less affected by emission in the off position. The CO temperatures 
are likely to be dominated by the gas heated by the PDR, and they therefore give only an upper limit to the 
mass averaged gas temperatures. The dust emission, however, is more optically thin and the submillimeter 
dust emission contains also emission from the colder dust inside the dense cloud ridge. This regime is more 



accurately measured by NH3, which shows kinetic temperatures of 25 K (Giisten & Marcaide 1986) along 
the cloud ridge. In the second entry for each offset lower gas temperatures are assumed, which are typically 
smaller than the CO temperatures. The two compact sources in the ridge have been substracted from the 
map in order to separate the dust cloud from the disk/envelope emission. However, the total mass values in 
Table 5 include also the contribution from LkHa 234 SMM 1 and LkHa 234 SMM 2. 



3.3. NGC 7129 FIRS 2 



NGC7129 FIRS 2 was first discovered by Bechis et al. ( 1978 ) as a cold FIR source associated with a 
^^CO column density peak south of LkHa 234. It was found to coincide with an H2O maser (Cesarsky et 
al. 1978 ; Rodriguez et al. 1980| ; Sandell & Olofsson 1981), suggesting that it is a young protostellar source, 



yet the source is invisible in the optical and near-IR. Additional FIR mapping by Harvey, Wilking & Joy 
(1984) confirmed the discovery by Bechis et al. ( |1978 ) and detected the source at both 50 and 100 /im, but 
not at 20 /im or shorter wavelengths. Edwards & Snell ( 1983| ) found that FIRS 2 was associated with a CO 
outflow, which further strengthened its association with a young, intermediate PMS star. Jenness, Padman 



& Scott (1995) included FIRS 2 in their survey of FIR cores in the vicinity of H2O masers and reported 



a flux density of 5.2 Jy at 800 /im and 30 Jy at 450 /im. Eiroa, Palacios & Casali ( 1998 ) proposed that 
FIRS 2 is an intermediate mass Class object, based on submillimeter photometry at 2 mm, 1.3 mm, 1.1 
mm and 800 /im as well as from high resolution reconstructed IRAS data. They detected FIRS 2 at 25, 60 
and 100 //m, but not at 12 /im. They fitted the SED of FIRS 2 with a standard isothermal graybody fit and 
found a /3-index of 0.9 and a dust temperature of 35 K corresponding to a total mass of 6M0. 

Our SCUBA data give a more accurate position for FIRS 2 (Table 3), which agrees well with the H2O 
position (uncertainty ± 10"). The source is resolved with a FWHM of < 3.7". The surrounding cloud 
core is also seen, but because most observations were made using relatively short chop throws ('^ 100"), 
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the emission of the cloud core is not weh mapped. The data is supplemented with maps and photometry 
made by us using UKT14 on JCMT. Our UKT14 flux densities are generally lower than those found by 



Jenness, Padman & Scott ( 1995 ) and Eiroa, Palacios and Casali ( 1998 ), but agree much better with the 
well-calibrated SCUBA data. We find flux densities of 0.48 ± 0.10 Jy, 1.18 ± 0.03 Jy, 1.69 ± 0.05 Jy and 
3.66 ± 0.22 Jy for 2 mm, 1.3 mm, 1.1 mm and 800 fim, respectively. The flux densities at 2 mm and 1.3 
mm are based on photometry and therefore also include a small contribution from the surrounding cloud. 
The 1.1 mm and 800 /im values are derived from Gaussian fits to calibrated maps, i.e. essentially the same 
way SCUBA data were analyzed. 

The isothermal fit (Fig. 5) predicts a /3 ^ 1 for dust temperature of 42 K. The fit underestimates 
slightly the FIR data points, which are observed with a large beam size and therefore also include emission 
from the surrounding cloud core. Our /3-index and mass (Table 4) agree rather well with the value derived 



by Eiroa, Palacios & Casali (199S) and demonstrate that FIRS 2 is a cold protostellar source with a 
luminosity of ~ 400 Lq . 



4. Discussion 

4.1. A Gas and Dust Comparison 

The ^'^CO emission (Fig. 1) and dust emission (Fig. 4) show broadly similar distributions. The 
molecular cavity is also a dust cavity, although the boundary appears broader in dust continuum emission 
than in ^■^CO emission. The presence of a dust cavity suggests that the interior is not filled with HI, but 
has been substantially cleared of dust and gas. Bechis et al. (1978) showed that radiation pressure from 
older stars in the cavity is capable of doing this. Moreover, recent 21cm HI observations by H. E. Matthews 
et al. (2001, in preparation) in the NGC 7129 region support this hypothesis, by showing that HI in the 
cavity has a low number density, of 85 cm~'^, and it is in pressure equilibrium with the molecular gas in the 
ridge (both having P/k 2x10^ cm^^K). 

In both ^^CO and dust, the mapped region is seen to contain two main clumps, one surrounding the 
star LkHa 234 and extending to the northwest, the second being about 2' to the south. Although the dust 
and gas maps are consistent on large scales, there are significant differences on smaller scale. The southern 
clump, centered on the peak NGC 7129 C04, has an extension to the east with no counterpart in the 
continuum maps. LkHa 234 SMM 1 (identified as NGC 7129 COl) dominates the emission in submillimeter. 
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while NGC 7129 C02 (SMM2 in Table 3) is almost equally bright in ^^CO and much fainter in dust 
continuum. The peak NGC 7129 C03 is not present in the dust emission maps. The two ^^CO peaks, NGC 
7129 C02 and C03, have essentially equal intensity, so the difference in continuum emission between the 
two is surprising. The answer may lie in the proximity of the star SVS 13 to SMM2. Radiation from SVS 
13 may be responsible for maintaining the dust temperature of ~ 30 K, which we find for this source. The 
projected distance of SMM2 from SVS 13 is about three times that of NGC 7129 C02 and therefore the 
stellar radiation field is less intense. However, the absence of dust emission at NGC 7129 C03 cannot be 
explained by temperature alone. It must have a lower density than NGC 7129 C02. 

A quantitative comparison of gas and dust is provided by the masses listed in Tables 2 and 5. The 
mass from the CO isotopomers was obtained from LTE expressions (Section 3.1), while the mass from the 



850 fim flux required values of /3, T^ and a gas/dust mass ratio of 100 (Hildebrand 1983) 



M{Mq) = 1.88 X 10"^ -kF X {-^TT^-f"^ X [e^^d - 1] x (3) 

where F is the flux (in Jy) at the wavelength A (in mm), d is the distance to the source (in kpc) and 
the mass opacity coefficient is 10 cm^/g at 250 /im. Two values for dust masses are given in Table 5, for 
different values of T and (3. We will refer below to those found with /3 = 2, which is the typical value 
expected from theoretical grain models (Hildebrand 1983| ). In view of all the assumptions that have been 



made, we consider the mass values derived from the CO and from the dust to be in good agreement. In 
the worst cases, the difference is a factor of three. In all these cases, the agreement could be improved by 
making defensible adjustments to the dust temperature. The results are also consistent with little or no 
CO depletion in the NGC 7129 ridge, perhaps due to the enhanced radiation field, which maintains grain 
temperatures above the 16 K CO freeze-out temperature. 



4.2. The nature of the submillimeter sources LkHa234 SMM 1 and SMM2 



Early studies of NGC 7129 (Bechis et al. |1978| ; Harvey, Wilking & Joy |1984| ) indicated that the FIR 
emission could be explained by dust heated by the B stars that illuminate the reflection nebulae. The FIR 
emission peaked on LkHa 234, which was also known to be a weak radio continuum source (Bertout & 



Thum 1982), associated with an H2O maser (Cesarsky et al. 1978 ) and driving a high velocity CO outflow 



(Edwards & Snell |1983[ ). The second FIR peak, FIRS 2 was also associated with an H2O maser and a high 
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velocity CO outflow, but lacked an optical counterpart (Section 3.3). Since LkHa 234 was associated with 
HH objects (Hartigan & Lada 1985| ), coincided with a free- free emission source at 15 and 22 GHz, with 



compact dust emission at 3 mm (Wilking, Mundy & Schwartz 198(: ) and drove a highly coUimated optical 



jet (Ray et al. 1990 ), this identification appeared rather secure. 

However, Skinner, Brown & Steward ( 1993| ) noticed that the free-free emission was offset by ^ 1'.'7 
from the optical position of LkHa 234 if they adopted the more accurate optical position given by Herbig & 
Bell (|198S| ) or by 2'.' 2 when compared with the astrometric position of Clements & Argyle ( |l984[ ). Near-IR 



polarimetric imaging by Weintraub, Kastner & Mahesh ( 1994 ) showed that there is a deeply embedded 
companion ~ 3" North- West of LkHa 234, which was invisible in total intensity maps at K band. High 



resolution mid-IR imaging by Cabrit et al. (1997) identified the embedded companion star, IRS 6, at 2'.'7 
NW of LkHa 234 and showed that it has a steeply rising spectrum, making IRS 6 as bright as LkHa 234 at 
17/im. IRS 6 therefore coincides with the radio free-free emission source, with two H2O maser spots, and 



with the 3mm dust emission (Wilking, Mundy & Schwartz 1986; Sandell, unpublished data). Cabrit et al. 



(1997) show that IRS 6 illuminates an arc-shaped refiection nebula with very red colors and that it drives 
an H2 jet at a position angle (p. a.) of 226°. The p. a. of the H2-jet differs from that of the optical [SII]-jet, 
which has a p. a. of 252° (Ray et al. |l990| ). Whether the much larger [SII] jet also originates from IRS 6 



is not clear (Cabrit et al. [19971) . Ahhough MitcheU & Matthews ( [I994D report that the p. a. of the CO 
outflow agrees with the optical jet, this is not true if the p. a. is determined from the symmetry axis of the 
large scale red CO outflow, which is ^ 230°. Neither does the p. a. of the optical jet agree with the location 
of the two HH objects CCD 32 and HH 103 SW of LkHa 234 and IRS 6, which lie in the blue-shifted lobe 
of the CO outflow. 

The submillimeter position for LkHa 234 SMM 1 determined from our 850 /im map (Table 3) falls 
between LkHa 234 and IRS 6, but since our positional accuracy is ^ 1", the emission could peak on 
either star. However, since the dust emission at 3 mm peaks on IRS 6, it is plausible to assume that the 
submillimeter emission also peaks on IRS 6. It therefore appears that it is the cold, heavily obscured IRS 6 
that drives the CO outflow and which is associated with the massive dust disk or envelope we see in the 
submillimeter. This disk-like structure is orthogonal to the near-IR H2 jet to within a few degrees. From 
our simple isothermal modeling (Section 3.2) we derive a luminosity of ~ 9 10^ Lq, which would make 
IRS 6 an early B star, i.e. a young Herbig Be star still heavily enshrouded by dust. It is clear that LkHa 
234 as weU as BD+65°1638 and even BD+65°1637 and SVS 13 contribute to the 50 and 100 ^m emission, 
but the submillimeter emission is totally dominated by IRS 6 and the second submillimeter source LkHa 
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234 SMM2. There is some dust emission in the molecular cavity as well, especially in the vicinity of 
BD+65°1638. This dust is hotter than the dust in the cloud ridge, which is also seen from the spectral 
index plot (Fig. 6). 

The submillimeter source LkHa 234 SMM2 does not have a known optical or IR counterpart. It is ~ 
12" to the east of SVS 13, which does not appear to be associated with any dust emission (Fig. 4). SMM2 



does not coincide with the CO hotspot seen by Mitchell & Matthews ( 1994 ) either, which is ^ 10" - 15" SE 
of the submillimeter source. SMM2 coincides within errors with the NGC 7129 C02 peak (Fig. 1) and is 



also apparent in the CS J=3— 2 map presented by Fuente et al. ( 1998 ) as a dense gas condensation. Two of 
our C^^O J=3— 2 spectra, at offsets 3 and 6 in Table 2, are close to SMM2, i.e. ~ 5" north and south of 
the submillimeter peak. They show slightly elevated excitation temperatures (see Table 2), but the higher 
kinetic temperature is more likely due to the FDR rim than to any embedded source. The linewidths, 
especially in C^^O, are narrow and there is no evidence for outflow. SMM2 therefore has the appearance 
of a pre-stellar core. The virial mass is ^ 8.2 M©, assuming a linewidth of 1 km s^^, and a diameter of ~ 
14" i.e. several times larger than the mass we derive from our submillimeter observations. However, it is 
clear that the the hot gas in the cavity has compressed the whole ridge in which SMM 2 is embedded. We 
also note that our derived value for /3 as well as our assumed dust temperature is rather uncertain and our 
mass estimate is at most accurate to a factor of two. SMM 2 could therefore either be a pre-stellar core, or 
a protostar. 



4.3. Dust properties 

To investigate dust properties across the region, a map of the spectral index a {Fn ^ v") was made, 
using the 850 and 450 /im deconvolved images smoothed to 14" resolution. Explicitly, 



a = log(FA,/^^Aj/log(A2/Ai), (4) 
where Ai and A2 are the two wavelengths used in our observations. 

We know that the maps are distorted by chopping into the cloud emission (Section 3.2). This will affect 
the 450 /im more than the 850 /im flux, because the error lobe is larger at 450 /im. The 450 /im map was 
therefore corrected by adding a constant level of 5% of the peak value to the map, before convolving it to a 
resolution of 14". The resulting map is shown in Figure 6. 
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The spectral index map has two trends. Firstly, a comparison with the 450 /im or 850 map shows a 
general anticorrelation between submillimeter continuum flux and spectral index. The highest flux positions 
are minima of a, and the value of a increases towards regions of weaker continuum emission. The a minima 
are seen at locations correlating with dense gas clumps, LkHa 234 SMMl, SMM2 and the flux peak near 
the offset (— 45'.'0,— 85'.'0). The curving ridge of 450 /zm emission (Fig. 4) can be followed throughout its 
length in the a map. Secondly, the map also shows a steep rise in a towards the molecular cavity. This 
large gradient in a occurs along the entire transition region. 

It is important to determine whether the observed increase in a is due to the change in temperature 
through the PDR, or in the dust opacity index, /?. Because F,^ oc Kv^Bi,{T4), the spectral index, a, and the 
dust opacity, (3, can be related as follows: 



where 7 can be thought of as a correction to the Rayleigh-Jeans form of the Planck function and is 
given by (Visser et al. 1998) 



An increase in a may bo due to an increase in the dust temperature, T, and/or in /3. The temperature 
is expected to be higher in the cavity because of the more intense radiation field. Our analysis in Section 
3.1 showed that the gas temperature is higher towards the ridge. Is it possible to understand the larger 
values of a only as due to higher dust temperatures? The answer is no because 7 approaches zero at large 
T so that the maximum value of a is /? + 2. Values of a larger than 4 require /3>2, even if T is large. While 
the expected increase in T^ toward the cavity will result in higher values of a, the increase in T^ must be 
accompanied by an increase in /3. 

Ideally, the above relations can yield the exact value of /3 if one knows the spectral index a and 
the temperature of the dust. In practice, the determination of the true value of /3 is restricted by the 

current uncertainties in the submillimeter fluxes, especially at short wavelengths, which will reflect into 
uncertainties in a. The dust will be more optically thick at 450 /im than at 850 /im (since the optical depth 
changes with frequency as r oc i'^). Thus at 450 /im we do not see as deeply into the densest hot cores as 
one does at 850 /zm. Another more severe restriction is our ability to accurately measure the submillimeter 



a; = /3 + 2 + 7. 



(5) 



7 = 1- log[(e^'^-<'/'=^ - l)/(e^-«o/'=^ - l)]/log(z.85o/z^45o). 



(6) 
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fluxes, especiafly at short wavelengths. As mentioned before, we tried to remedy this situation by adding a 
constant zero-level to the 450 fim flux, but this is only a flrst order correction. In reality, the dust structure 
will not be flat, but will have its own substructure, similar to that of the molecular cloud. Therefore, the 
numerical values of a in regions of low flux are highly uncertain. 

In summary, we have confidence in the values of a ~ 3 that we find at high flux levels. We also believe 
in the trend of increasing a towards lower flux values. In the dense cores we flnd /3 < 2, as suggested by 
the minima in the spectral index plot as well as our dust flts in Figure 5. Since small grains produce a high 
/3-index and large grains lead to smaller values of /3, our results are consistent with the growth of grains in 
regions of high density. 

These present observations of NGC 7129 support the conclusions of other similar studies that high 



density clumps or cores have smaller values of /? than does the surrounding gas (Visser et al. 1998; 
Johnstone & Bally 1999 ; Sandell et al. 1999| ). Values of /3 ^ 2 outside the flux peaks have been also 



reported in reflection nebulae surrounding Herbig Ae/Be stars (Whitcomb et al. 1981 ). Our a map goes 
further, however, in showing a general correspondence between flux and a over an extended area, and in 
showing a steep increase in a (and, therefore, in /3) along an extended PDR. 



5. SUMMARY 

New high resolution maps of the NGC 7129 ridge have been obtained, both in submillimcter continuum 
emission and in ^■^CO J=3— 2 line emission. The main results of this study are: 

1. A sharp boundary between the molecular ridge and cavity is clearly seen in the ^^CO map. A 
previously detected narrow fllament of optical and near infrared emission follows closely the cloud-cavity 
boundary. Strong emission in the v=l— S(l) line of H2 from the fllament, if due to fluorescence, supports 
the interpretation of the boundary as a PDR. 

2. The ^^CO emission consists of two main cloud cores connected by a bridge of gas. The northern 
cloud core has three ^'^CO peaks, while the southern one has a single extended peak. 

3. The 850 /im and 450 /im maps show a very good morphological agreement with the ^'^CO map. 
Both data sets show a sharp boundary (cloud rim) toward the cavity where the two young stars BD +65 
1637 and BD +65 1638 are located. 
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4. Three compact sources are seen in the continuum maps (SMMl, SMM2, and FIRS 2 in Table 3). 
Isothermal graybody fits to the measured spectral energy distributions have been made, by combining our 
SCUBA data with earlier observations. The resulting masses, dust temperatures, luminosities and dust 
emissivities are given in Table 4. The SMM 1 continuum peak is identified with the deeply embedded 
infrared source IRS 6, which according to these present fits is a young Her big Be star. 

5. For nine offsets around the LkHa 234 star, the ^^CO/C^*0 intensity ratio was used to obtain the 
mass of gas in an 14" beam. The masses derived from CO are found to be consistent with masses derived 
from the dust continuum emission. Within our analysis uncertainties the results are consistent with little 
or no CO depletion in this region. 

6. A map of the 850/im/450/im spectral index, a, shows a general correlation between the a minima 

and the submillimeter continuum flux peaks. The low values of a in the peaks also correlate with low values 
of dust opacity indices, /3, in the range 1 to 1.6, derived for the compact source, suggesting grain growth. 
The spectral index a rises steeply towards the molecular cavity, along the entire length of the boundary 
region. Consequently, the dust opacity index in this region is also high, (3^2. This suggests very hot, 
small size grains in the PDR ridge, whose mantles have been evaporated by the intense UV radiation. 

The James Clerk Maxwell Telescope is operated on a joint basis between the United Kingdom Particle 
Physics and Astronomy Research Council (PPARC), the Netherlands Organization for the Advancement 
of Pure Research (ZWO), the Canadian National Research Council (NRC), and the University of Hawaii 
(UH). This research was partly supported by an operating grant from the Natural Sciences and Engineering 
Research Council of Canada. 
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Table 1. Physical parameters of the ^^CO peaks. The peak antenna temperature, the line width and line 
velocity center were determined with a Gaussian model best fit of the spectral lines. 



Source 


a (2000) 


5 (2000) 


T*A 


Av 


^center 




^/i m 5^ 


(O , 


(K) 


(km s-i) 


(km s-i) 


NGC 7129 COl 


21 43 06.7 


+66 06 57.0 


13.08 


2.96 


-10.30 


NGC 7129 002 


21 43 03.5 


+66 07 10.0 


15.02 


1.21 


-10.16 


NGO 7129 003 


21 43 02.0 


+66 07 28.0 


15.89 


1.24 


-10.12 


NGO 7129 004 


21 42 58.5 


+66 05 30.0 


17.55 


1.44 


-10.39 
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Table 2. Measured properties of the gas for the nine offsets where both C^^O and ^^CO observations are 

available? 



No. Offset T*a{^^CO) T*a {&^0) ns Mg„, 

(arcsec) (K) (K) (K) ( M0/I4" beam) 



1 


(52.5, 37.5) 


3.84 


1.32 


4.13 


13.1 


3.6 


2 


(-5.0, 13.0) 


12.57 


4.43 


4.27 


29.2 


7.3 


3 


(-20.0, 25.0) 


15.73 


7.68 


6.69 


34.7 


12.8 


4 


(5.0, -8.0) 


13.13 


3.87 


3.35 


30.7 


11 


5 


(0.0, 0.0) 


12.8 


4.18 


3.85 


29.8 


10.7 


6 


(-22.5, 10.0) 


14.97 


4.96 


3.92 


33.6 


7.2 


7 


(-45.0, -85.0) 


17.54 


5.70 


3.82 


38.3 


8.3 


8 


(-30.0, -100.0) 


13.04 


3.25 


2.63 


31.5 


6.7 


9 


(-20.0, -110.0) 


7.84 


2.97 


4.71 


20.6 


6.4 


*The 


last two columns 


in the table 


were calculated assuming 


LTE 


conditions 


and that the 



thermalized. 



Table 3. Positions and deconvolved Gaussian sizes of compact submillimeter sources. 



Source a(2000.0) 5(2000.0) Ga x Ob P.A. 

C' ™ ^) (° ' ") (arcsec) (degree) 



LkHQ!234 SMMl 21 43 06.76 +66 06 56.0 6.1 x 4.9 -46 

LkHQ!234 SMM2 21 43 03.20 +66 07 13.1 14.6 x 9.6 -67 

NGC7129FIRS 2 21 43 01.51 +66 03 24.2 3.9 x 3.7 n.a. 
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Table 4. Physical properties of the deconvolved compact submillimeter sources. 



Source 


S850 


S450 







Mtot 






(Jy) 


(Jy) 


(K) 




(Mo) 


(L0) 


LkHa 234 SMM 1 


3.12 


20.7 


45 


1.24 


7.6 


8.8 X 10^ 


LkHa 234 SMM 2 


0.73 


6.2 


30 


1.62 


4.8 


1.4 X 10^ 


NGC 7129 FIRS 2 


3.35 


18.1 


42 


1.04 


7.0 


4.3 X 10^ 
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Table 5. Dust properties for the nine C^'^O positions. The second entry for each offset uses the cloud 
emission only (when applicable) and a value of 2 for the /3 index. 



Offsets SssO/nm S450^m Trf /3 Mtot 

(arcsec) (Jy/14" beam) (Jy/14" beam) (K) (M0/I4" beam) 



(52.5, 37.5) 


0.26 


1.39 


13.1 


1.71 


6.5 




0.26 




13.1 


2.0 


9.3 


(-5.0, 13.0) 


0.74 


6.36 


29.2 


1.70 


5.5 




0.43"* 




25.0 


2.0 


8.2 


(-20.0, 25.0) 


0.45 


5.18 


34.7 


2.06 


4.1 




0.26'" 




25.0 


2.0 


4.6 


(5.0, -8.0) 


1.44 


11.6 


30.7 


1.58 


8.6 




0.47^" 




25.0 


2.0 


8.5 


(0.0, 0.0) 


3.05 


23.1 


29.8 


1.50 


13.9 




0.35^" 




25 


2.0 


11.2 


(-22.5, 10.0) 


0.47 


4.15 


33.6 


1.68 


2.8 




0.19"* 




25.0 


2.0 


4.3 


(-45.0, -85.0) 


0.42 


3.33 


38.3 


1.46 


1.6 




0.42 




25 


2.0 


5.9 


(-30.0, -100.0) 


0.19 


2.24 


31.5 


2.14 


2.2 




0.19 




25.0 


2.0 


2.5 


(-20.0, -110.0) 


0.13 


1.83 


20.6 


2.67 


7.4 




0.13 




20.0 


2.0 


2.3 



''Cloud emission only, emission from LkHa 234 SMM 1 and LkHa 234 SMM 2 has been subtracted. 
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Fig. 1. — -'^^CO 3=3—2 integrated intensity map superimposed on an K' image of the region (Hodapp, 1994). 
The map was centered near the star LkHa 234, with the center coordinates a(2000) = 21'' 43™ 06". 08 and 
5(2000) = 66° 06' 56'.'29, and extending from -100" to +100" in a offsets and from -150" to 75" in S 
offsets. The crosses indicate the nine offsets with complementary C^^O J=3— 2 data which have been chosen 
for analysis (the numbers correspond to Table 2 entries). 

Fig. 2. — Composite 850 fim map of LkHa 234 and NGC7129 FIRS 2 in grey scale overlaid with contours. 
The 850 /im flux density is contoured from 50 mJy/beam in logarithmic steps equal to a factor of 10*^'^ of 
the previous level. The size of the beam, 14"5, is indicated at the bottom right of the plot. Near-IR sources 
are labeled by stars and H2O masers by triangles. 

Fig. 3. — ^'^CO and C^^O J=3— 2 spectral lines shown superimposed for the nine offsets marked in Figure 1 
and listed in Table 2. The thin lines correspond to ^^CO emission and the thick lines to C^*0 emission. 

Fig. 4. — Contour maps of deconvolved images at 850 /im (14" beam) and 450 /im (8" beam). For 850 /im 
the contour levels start from 100 mJy/beam and go in steps equal to a factor of lO'''^ of the previous level. 
For 450 /xm the contour levels start at 250 mJy/beam and go in steps of 10"-^^ of the previous level. 

Fig. 5.— Dust fits for LkHa 234 SMM 1 and NGC7129 FIRS2, the latter data being offset by a factor 
of 100. The data points include our SCUBA observations, previously unpublished UKT14 data (Sandell, 
private communication), as well as IRAS and KAO data (Bechis et al. 1978; Harvey, Wilking & Joy |l984 ). 



The latter observations were not considered in the fit, since they are likely to include hot dust from the 
surrounding reflection nebulosity, because of their large beam sizes. Fitted data points are marked with 
filled triangles, all other by filled squares. 

Fig. 6. — Spectral index a map [F^ ~ v"^) obtained from the (14" beamsize) 850 ^m and 450 ^m maps. 
Greyscale contours range from 3 to 5, with a step of 0.25, where dark color implies high spectral index. The 
optical positions are marked with stars and the LkHQ;234 SMM 2 source with a triangle. 
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